Ground state of a tightly bound composite dimer immersed in a Fermi Sea 
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In this paper we present a theoretical investigation for the ground state of an impurity immersed 
in a Fermi sea. The molecular regime is considered where a two-body bound state between the 
impurity and one of the fermions is formed. Both interaction and exchange of the bound fermion 
take place between the dimer and the Fermi sea. We develop a formalism based on a two channel 
model allowing us to expand systematically the ground state energy of this immersed dimer with the 
scattering length a. Working up to order a"^, associated to the creation of two particle- hole pairs, 
reveals the first signature of the composite nature of the bosonic dimer. Finally, a complementary 
variational study provides an accurate estimate of the dimer energy even at large scattering length. 

PACS numbers: 03.75.Ss; 05.30.Fk; 34.50.Cx 
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I. INTRODUCTION 

As demonstrated by Bardeen-Cooper and SchriefFer, 
superconductivity arises from the pairing of electrons 
with opposite spins into Cooper pairs A natural 
extension of their work to the case where the two spin 
populations are imbalanced was proposed by Clogston 
and Chandrasekhar 0, Q : they suggested that when a 
magnetic field is a applied, the existence of a pairing gap 
in the electron energy spectrum could prevent spin polar- 
ization as long as the Zeeman shift was smaller than the 
gap, a threshold known as the Clogston-Chandraskhar 
(CC) limit. It was later suggested that the superfluid 
state could survive beyond the CC limit, into the form of 
a non-homogeneous superconducting state known as the 
FFLO (Fulde Ferrell Larkin Ovchinnikov) state [1, 
However, most superconductors behave like nearly ideal 
diamagnetic compounds (Meissner Effect), which forbids 
magnetically induced spin polarization in the bulk, see 
Ref. 6] for a short review. As a consequence, theses the- 
ories were investigated experimentally only very recently 
in a series of works performed at Rice [3] and MIT Q with 
ultra-cold Fermi gases trapped in optical potentials. Al- 
though some debate on the structure of the normal com- 
ponent persists between the two groups, both observe a 
shell structure in the density profile, with at center a fully 
paired region consistent with the CC scenario of robust 
fermionic superfluidity. 

The main discrepancy between the two experiments 
lies in the polarization of the normal component: In- 
deed, while Rice's group observed that the outer rim was 
composed exclusively of majority atoms, the normal com- 
ponent obtained at MIT was only partially polarized and 
contained also particles of the minority spin species. Re- 
cent theoretical work demonstrated that this latter obser- 
vation was compatible with the homogeneous phase dia- 
gram of a strongly interacting Fermi gas if one assumes 
the validity of the Local Density Approximation . In 
particular, it was shown that some of the most salient 
features could be understood fairly accurately from the 



study of the simpler problem of an impurity immersed 
in a Fermi sea [13, ]1M ■ unitarity (scattering length 
a — oo), it was demonstrated that the impurity could 
be described as a quasi-particle dressed by particle hole 
excitations of the background Fermi sea [ill . [T^ . [Tsl . [l^ . 
Comparison between simplified variational models and 
Monte-Carlo simulations have in addition shown that 
even in this regime of strong correlation, a single particle- 
hole excitation was sufficient to capture quantitati vely 
the properties of this so-called Fermi polaron [13, Il5j |. 
However, the Fermi-polaron picture is valid only for a 
negative, and close to unitarity 1/A:i?|a| ^ 1. In the 
l/kpa < 1 regime, Monte-Carlo simulations have re- 
vealed that the Fermi-polaron was not the ground state 
of the system anymore [l^. Indeed, for kpa < 1.11, the 
ground-state is now described by a bosonic dimer inter- 
acting with the background Fermi sea where the atom- 
dimer scattering length [l^ Udm — 1.1786 a was first cal- 
culated by Skorniakov and Ter Martirosian [l3|. 

In this paper, we use a combination of perturbative 
expansion and variational calculation to draw a simple 
and intuitive picture of the molecular sector of the impu- 
rity problem and extend diagrammatic calculations pre- 
sented in [H, [l^: we show that for a > 0, the molecular 
impurity shares several features with the Fermi-polaron 
describing the impurity for a < 0. In particular, we 
show that the creation of a single particle-hole pair in the 
Fermi sea provides an accurate quantitative description 
of the system. The paper is organized as follows: Sec. |TT] 
presents the main results of the paper. More technical 
aspects of this work are then relegated to following sec- 
tions, Sec.[in]for the perturbative calculation and Sec.HVl 
for the variational approach. Sec. IVl concludes. 



II. MAIN RESULTS 

Let us consider an ensemble of spin F — 1/2 fermions 
(for instance ^Li in its hyperfine ground state), where 
all particles but one are polarized in mp = —1/2. At 
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FIG. 1: Our variational scheme takes into account only a 
restricted number of particle-hole excitations. First row: no 
particle hole excitation. Second row; one particle-hole excita- 
tion of momenta q < kp and k > kp. Third row: two particle- 
hole excitations of momenta q, q' < kp and k, k' > kp. 



low temperature, short range interatomic interactions be- 
tween identical fermions are suppressed by Pauli exclu- 
sion principle, and the system can therefore be described 
as a single impurity interacting with a Fermi sea of non- 
interacting fermions. We note a the scattering length 
between particles of opposite spins. The case a < and 
the vicinity of the Feshbach resonance (a oo) have 
been described in previous works on the Fermi-polaron 
[O, [3 ' ^^'^ here we concentrate on the molecular 
sector corresponding to a small and positive. In this par- 
ticular regime the picture is relatively simple, since in the 
absence of a surrounding Fermi sea, the two body poten- 
tial possesses a deeply bound state of size ~ a and energy 
E = —h^/ma^. When one adds a spin mp = +1/2 atom 
to the Fermi sea, it will pair up with a majority atom 
to form a deeply bound dimer. In the regime kpa <C 1, 
the size of this molecule is much smaller than the inter- 
particle distance, and we can describe the dimer as a 
point like boson of mass 2m. In this picture, the energy 
shift associated with the addition of the impurity is then 
at leading order in kpo. 



where the second term corresponds to the removal of 



one majority fermion from the Fermi sea to form the 
bosonic dimer. The next order in the kpa expansion 
comes from the interactions between the bosonic dimer 
and the surrounding Fermi sea. Indeed, if the point-like 
boson picture is correct, then one should expect a mean- 
field energy shift AE = ga.dn, where n is the density 
of majority atoms and ^ad is the atom-dimer coupling 
constant given by gad = 37^?^^aad/n^ associated with the 
atom-dimer scattering length Oad ~ 1.1786a [Hi, Al- 
though this scenario has been confirmed by Monte-Carlo 
simulations [T5| . the analytic calculation of the mean 
field shift is not trivial for two main reasons. First, the 
atom-dimer scattering length is obtained by solving the 
three body-problem, but here antisymmetrization of the 
global wave-function correlates automatically the major- 
ity atom bound in the dimer with the surrounding Fermi- 
sea, and turns the calculation from three to many-body. 
Second, mean field contributions are usually obtained us- 
ing a perturbation expansion based on the existence of 
a small parameter (here kpa), and in our situation, this 
scheme is due to fail. Indeed, in absence of interactions 
(a = 0), the system is just an ensemble of non inter- 
acting fermions, that cannot form any dimer, which is 
contradictory with the intuitive picture we drew above. 
To circumvent this latter point, we decided to work in 
the two-channel picture [lo, HH , where the short-range 
two-body bound state responsible for the Feshbach res- 
onance is explicitly included in the model as a bosonic 
degree of freedom, of mass 2m and bare binding energy 
Eq. In second quantized form, the Hamiltonian reads 

H= efe4,a«fc,. + II(£^o + e/f/2)5W 

fe.<T=T,i K 

+ X! (yK"-k+K/2,]'^-k+K/2,l + ^■^■) > 

k,K ^ ^ 

where is a quantization volume and is the matrix el- 
ement coupling fermionic and bosonic degrees of freedom. 
If ri, is the typical size of the bare molecule, then the 
width of Afc is ~ l/rb. The fermionic operator ak,a de- 
scribes a momentum-fc atom with spin (hyperfine state) 
a. Atoms in the two different spin states, a =t,i, are 
coupled via a molecular closed channel state described by 
the boson operator with bare energy at rest Eq and 
momentum K . With an even (and very weak, see below) 
/c-dependence for the coupling strength A^, the resulting 
atomic interaction has an s-wave character. In the uni- 
versal regime where the properties of the system depend 
only on the scattering length, no other ingredient is nec- 
essary, as long as Eq and Aj, are chosen to reproduce the 
actual interatomic scattering length. The Pauli principle 
imposes vanishing s-wave interaction for fermions in the 
same spin channel and do not require coupling between 
same spin particles. 

Even in the case where a single spin down particle is 
immersed in a Fermi sea of spin up particles, this Hamil- 
tonian does not have a complete analytical solution. One 
needs to rely on approximations to get some insights on 
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the behavior of the system. It can be assumed, in analogy 
with the Fermi-polaron problem, that the most salient 
features are obtained by taking into account the exci- 
tations of a small number of particle-hole pairs of the 
background Fermi sea. In the work presented here, we 
take into account up to two particle- hole pairs, which re- 
veals effects due to the composite nature of the bosonic 
dimer. The variational state that we have studied thus 
takes the form (see also Fig. [1]) 
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where means that the sum over the majority parti- 
cle (resp. hole) wavevectors k (resp. q) is restricted to 
k > kp (resp. q < kp), and \FS) is the noninteracting 
Fermi sea of majority atoms, in the absence of minor- 
ity particle or short range molecule. Solving exactly the 
equations in this restricted subspace is a numerical dif- 
ficult task outside the scope of this work. We therefore 
simplify the problem in two directions: one leads us to a 
systematic perturbative expansion in kpa, the other one 
to a variational prediction for the molecule energy. 

We first expand the variational equations in power of 
kpa. In fact, as it was stressed in (T3| . an expansion 
in kpa is equivalent to an expansion in the number of 
excited particle-hole pairs. In particular, the variational 
space ^ is sufficient to derive the exact and systematic 
expansion of the immersed molecule energy up to order 
(kpa)^. To be more precise, extending the variational 
space by allowing additional electron-hole pairs would 
lead to corrections that are higher orders in kpa. 

Interestingly, the mathematical structure of the per- 
turbative equations exhibits a separation of length scales, 
kpa <^ 1, which can be given a simple physical interpre- 
tation. (pk,q is seen as the atom-dimer wavefunction when 
the Fermi sea is composed of a single majority fermion. 
It presents a singularity ^ \/k^ at small k corresponding 
to the general large distance decay ^ 1/r for the scat- 
tered wave of a short-range potential. The prefactor gives 
the renormalized atom-dimer scattering length. Due to 
this singular behaviour, the wavefunction ^pk.q is modi- 
fied separately at the two length scales a and Xp = 1/kp 



in the presence of the Fermi sea. It is first modified at 
small wavevectors ~ kp, i.e. close to the singularity, the 
stronger effect being Pauli blocking which restricts al- 
lowed wavevectors k to lie outside the Fermi sea. The 
corresponding contributions to the molecule energy are 
the same as in the point-like boson case, like the first two 
corrections in Eq. ^ for example. They involve large dis- 
tances ~ Af 3> a over which the composite nature of the 
molecule is not visible. In contrast to that, corrections to 
Vk,q for larger wavevectors ~ 1/a reveal the composite 
structure of the molecule and change the value of as in 
Eq. ©. 

Using the simplifications permitted by the length-scale 
separation mentioned above and working out the expan- 
sion in Xp = kpa to third order yi^i^'^^ following ex- 
pression for the molecule energy [2g| (see Eq. ((55)) ') 
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(In 2- 3/8) 



(3) 



where Oad — aad/o- ^ 1.1786 is the ratio between the 
atom-dimer and atom-atom scattering lengths, and ~ 
0.0637 is a numerical coefficient obtained from the for- 
malism described in Sec. IIIIl As mentioned above, the 
first term of the expansion corresponds to the binding 
energy of a single molecule. The next two terms are 
obtained by taking into account single particle-hole ex- 
citations and depend only on the atom-dimer scattering 
length. They are identical to the energy shift obtained for 
a point-like boson immersed in a Fermi sea (2^ and the 
composite nature of the dimers is only revealed by the 
last term of the expansion, which is calculated by tak- 
ing into account two particle-hole excitations. We stress 
again that, although the starting point is a variational 
form ([2]), Eq. ([3]) is the exact low density expansion of 
the molecule energy. For comparison, let us mention that 
the energy obtained for a point-like boson immersed in a 
Fermi sea with boson-fermion scattering length a^d also 
expands as Eq. ([3]) - from which the molecule energy in 
vacuum has been subtracted - but with ~ 0.00025 ai^^j 
instead of as ~ 0.0637 in the composite boson case. 

It has been pointed out in [l^ that the exact energy 
of the molecular impurity was actually very close from 
the mean- field correction, up to kpa ^ 1. Due to its 
perturbative nature valid only for kpa <C 1, the anal- 
ysis presented above cannot explain this intriguing fea- 
ture. To address this particular issue and acquire some 
insight on the strongly interacting regime where the per- 
turbative expansion diverges, we develop a variational 
treatment of the impurity problem. In order to get a 
simple and tractable calculation, the general variational 
form ^ is simplified by allowing only one particle-hole 
pair excitation and by suppressing the q dependence in 
the wavef unctions. Namely, the variational state takes 
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FIG. 2: Perturbative expansion of the energy of a molecular 
impurity immersed in a Fermi sea. Square: Monte-Carlo re- 
sult of ifSj]. Black solid line: Free molecule E = —h^/ma^. 
Red dashed line: single particle hole excitations correspond- 
ing to the expansion up to a^. Dotted blue line: expansion up 
to taking into account a second particle-hole excitation. 
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IV-) = (/36it=o + E'^^"li"-fe.T 
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FIG. 3: Variational calculation of the energy of a molecule 
immersed in a Fermi sea. Solid line: Variational calcula- 
tion including formation of a single particle hole pair. Points: 
Monte-Carlo calculation of [Tsj . 

by step by increasing the variational space size up to two 
particle-hole excitations. The two-body problem is con- 
sidered in nil Al which helps us to connect the parameters 
of our starting Hamiltonian ([T|) to the scattering length a 
in the single-channel limit. A trivial ansatz extending the 
two-body wavefunction is proposed, which exhibits the 
Pauli blocking effect. Elaborating on this ansatz, single 
particle-hole excitations are added in IIII Bl yielding the 
first (mean-field type) and second order corrections to 
the molecule energy. The third order correction requires 
two particle-hole excitations that are considered in IIII CI 
and allows us to recover Eq. ([3]). 



As shown later, this particular ansatz gives the correct 
expansion up to mean-field term. In addition, it can be 
extended to strong interactions, and the energy remains 
finite even at unitarity. Working out the variational equa- 
tions, we reduce the problem to a set of two integral 
equations that are solved numerically, see Sec. IIVI The 
result is displayed in Fig. [31 We observe a remarkable 
agreement between our simplified integral equations and 
the Monte-Carlo simulations. It demonstrates that, just 
like the Fermi-polaron problem, single particle-hole ex- 
citations provide an accurate description of the physical 
properties of the system. 



III. PERTURBATIVE EXPANSION 

We shall detail in this Section our perturbative calcu- 
lation for the immersed molecule energy. A consistent ex- 
pansion at a given order 'mkpa requires a corresponding 
minimal number of particle-hole excitations in the vari- 
ational space. We will therefore do the calculation step 



A. Two-body properties 

Insight on the model ([T]) can be gained by solving 
the two-body problem. We look for a bound state of 
energy E — —^'-^ with the ansatz \^) = (/?fo^^o + 
J^k^^a^e l^^-k 1)1^)- Schrodinger equation {H — 

E)\^p) =0 gives simple coupled equations 

(g2 + A2)/?+^ AfeAfc = 0, Akf3+{k'' + X^)Ak = 0, (5) 

k 

with Afe — { j^2"yY j and Eq = . The second equa- 
tion allows to express Ak as a function of f3. Replacing 
Ak in the first equation yields an equation for the bound 

state 

^o + ^'-EpTa^ = 0- (6) 

k 

Although Afe has a very weak k dependence, the integral 
in Eq. ^ has an ultraviolet divergence if this dependence 
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is crudely neglected. We therefore define X]fc w ~ "1^' 
then add and substract this integral to Eq. Afc is 
replaced by Aq for the resulting convergent integral. A is 
finally solution of a second order equation 



a 



where we define 

f 

a 



1 Airh^En 



niAl 



0, 



R* = 



47rr 



(7) 



(8) 



where g 



The first term in Eq. ^ 



the 



molecule energy in vacuum, is the correct leading result 
in a small kpi expansion. The second term, associated 
with Pauli blocking, corresponds to a mean field term 
treated within the Born approximation and may naively 
be interpreted as a first sign of the composite nature of 
the dimer. However, as we will see, this result is not 
exact, and the correct value of the mean field term re- 
quires the resolution of the three body problem. This 
can be improved by including particle-hole excitations in 
the variational wavefunction as we will show next. 



i acts here as a short-distance cutoff length, it is on 
the order of the van der Waals potential size. The two- 
channel model Eq. ([T]) is therefore only valid if i?* ^ £ 
and a ^ £ for describing the two-body bound state. The 
two-body scattering problem can also be solved using the 
same ansatz. One concludes from its solution that a is 
the two-body scattering length and R* = —re/2 where 
Te is the effective range. The universal regime, or one- 
channel zero range model, is recovered with i?* — > that 
is Afe ^ oo at fixed value of a. In particular, Eq. ^ 
gives then A = f /a as expected. 

It is instructive to first study a simple ansatz where 
the dimer state and the spin J, Fermi sea are not coupled, 
namely 



(9) 



B. First order 

Applying two times the Hamiltonian ([T]) on the simple 
previous ansatz Eq. ([5]), one sees that the space spanned 
by the variational wavefunction can be extended to in- 
clude states with a molecule and a single particle-hole ex- 
citation from the Fermi sea. States with a closed channel 
boson and a single particle-hole excitation should also be 
included as an intermediate step. The improved ansatz 
for the wavefunction reads 



\ fe 



k.q 



(13) 



\FS) denotes the | Fermi sea with Fermi wavevector 
kp and density n = kp/Gn^ [13]. The Pauli principle 
imposes that wavevectors below kp are blocked in the 
FS and do not participate to the molecular state. This 
ansatz is certainly not an eigenstate of the Hamiltonian 
Eq. The calculation that follows is therefore vari- 
ational in essence. Projecting the Schrodinger equation 
{H—E)\'i(j) = onto the two relevant sectors, two coupled 
equations are derived 



k,k' ,q / 

where, thanks to the linearity of Schrodingcr's equation, 
we can take (3 = \. Apart from the additional q depen- 
dence, Eq. (|13p resembles the variational state Q. The 
variational space has four different sectors. Projecting 
the Schrodinger equation on the first two sectors, we get 



\k\>kF 



0, KkP+{k^ + \^)Ak = 0. 

(10) 

Apart from the wavevector summation restriction, they 
are identical to Eqs. ([5]). The spin | Fermi sea enters 
only through Pauli blocking at this stage. The energy 
E is measured from the Fermi sea energy. Solving the 
system in the universal limit R* — > 0, we find 



^ _ 1 2kp 
a TT 



1 ^ . 



(11) 



First, we can consider the dilute limit kpa <ti 1, which 
yields 



E 



+ 2gn + 0{gn{kpaf), 



(12) 



(9o + A')+ Y AfeAfc = 

\k\>kF 

Afe + (fc2 + A2)Afc + ^ 

\q\<kF 



(14a) 

k^k,q = 0. (14b) 



The aim of this calculation is to determine how the 2gn 
term (mean field) in Eq. (|12p is modified by the improved 
ansatz Eq. (|13p. Hence the q-dependence of (pk^q is not 
relevant because it leads to corrections of higher order in 
kpa. For R* ^ and kpa ^ 1 the modified bound state 
equation is found to be 



1 2k'j-^CL 2/ct^ V A 
A= ^ + 



37r 



37r -^^ + 1/02' 

k 



(15) 



where cph still has to be determined. 
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Again we neglect the qr-dependence when projecting 
on the last two sectors. Two coupled equations are de- 
rived where the anticommuting properties of fermionic 
operators are essential 



2^ ^k'+k/2^k.,k' 



ql + \^ + -k^ ) fk = -AoAfc = 



^^0 



(16a) 



P + A2 



Afc'+fc/2 {Vk - fk') + 2 (A^ + fc2 + _^ fc . fc') = 0. 

(16b) 

With an additional source term on the right side of 
Eq. p6ap . this system of equations is the same as the 
one derived in the three-body problem [2lj. The second 
equation is solved readily in '^k,k' and inserting this re- 
sult in the first equation, we obtain for i?* ^ 



47r 



ifk 



fc2 + A2 ' 

with the usual three-body kernel [l^ [l3| 



1 



+ 



27r2 



Vk' 



k'^ + k^ + k'-k + A2 



At this order, it is consistent to take A = 1/a. We proceed 
further and make contact to the atom-dimer scattering 
problem. Writing ipk = ^^ao(fca) yields the integral 
equation 



1 



au — m 

TT In uu' 



'0 

3 



1 



,'2 



ao(u) 



1 + 1*2 + u'2 — uu' 
1 



(19) 



4 1 + ^1 + 3u2/4 1 + ' 



identical to the one that determines the atom-dimer scat- 
tering length [H, In particular ao(0) = ^ ~ 
1.1786... 

The last term in Eq. is finally obtained by taking 
the fc ^ Umit in Eq. (fT7)l . 



E 



•i^k 



fc2 + l/a2 



1- o«o(0) 



(20) 



This contribution is positive so that the molecule energy 
is lowered by improving on the simple ansatz. This is in 
fact a consequence of a general variational principle which 
states that the ground state energy can only decrease 
when the parameter space is extended. Inserting Eq. ([^D]) 
into Eq. ([T5)) leads to a final expression for the low density 
molecule energy 



E 



(21) 



where the atom-dimer coupling constant has been defined 
as 



5ad 



2TTh Qad 

(2m/3) ' 



(22) 



Interestingly, the expansion can be pushed to the next 
order within the same variational space. Hence, the q 
dependence of ipk,q becomes relevant. Skipping details, 
the result for A reads 



|D|<a;F \u\>xp I 

where xp = kpa and stands for J .^Zi . (pu^v denotes 



a rescaled atom-dimer wavefunction 
solution of the integral equation 



(24) 



1 



(17) TT^ 



27r2 J\u'\>XF 1 + '"'■^ + ^2 + u'u — v{u' + u) 




3^2 U ■ V \ 1 

1 H : : 1 'Pu,v = 



(25) 



4 4 2 



1 + m2 



(18) For = 0, we recover Eq. with (fu,v — aQ{u)/i 



and therefore Eq. (PT|) for the energy. A first idea to 
treat Eq. (|25|) is to expand directly in powers of v and 
solve it order by order. Note that Eq. (^5)1 involves an 
average over the direction of v. This would imply that 
the first non- vanishing correction goes as v'^ and therefore 
Xp for the energy. However, this approach fails because 
the solution of Eq. (^5]) is singular as u — + 0, if>u^v — 
ao(0)/u2 for t> = 0. 

A correct treatment requires a proper description of 
(fu,v as a function of v for u ^ 0. We therefore write 



asi{u,v) 



m2 — ^2^3 _ 2u, • i>/3 ' 



(26) 



which translates Eq. (|25|) to an integral equation for 
asi{u,v). The second term in the l.h.s. of Eq. ([^5)1 can 
then be expanded in v. For the first term in the l.h.s. of 
Eq. (|25p . we add and substract the kernel for v = and 
Xp = 0. The resulting integral corrections are simplified 
by changing variables, u = vx and u = xpx, followed 
by an expansion in v and xp. Keeping only first order 
corrections, the integral equation reads 

du' KQ{u,u')asi{u' ,v) = 

(27) 

l + asi{0,v) [^^~^Mv/^f)] 
1 + u^ 

with the same kernel Kq{u,u') as Eq. p^ . We have 
defined 



Ji{s) = y j 



l + 2x/y 



yy3-2xy/3- 



(28) 
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The solution of Eq. ([271) is then simply 



asiiu,v) 



(29) 



and can be expanded to first order in xp and v. The 
first outcome of this calculation is that the singular be- 
havior of ifiu.v at small u indeed modifies the xp (and v) 
expansion. 

The solution ((29| can be included in Eq. ((23)l . Again 
we add and substract the solution at vanishing x p and v 
and rescale integral corrections with v and x p . The final 
result for the energy is 



1 



4«o(0) _ 2x%al{0) 
2tt 7r2 



(In 2- 3/8) 



which gives the next order xp = kpa correction to 
Eq. ([2T|l . We have used the following result 



/■I 7 4 

/ dss3ji(s) = - - -ln2. 
Jq 6 3 



(31) 



A few conclusion can be drawn from this result ([50)1 . As 
we have mentioned already, the first correction to the 
mean field result (PT|) is linear in xp and not quadratic. 

t 2 

In fact, subtracting the dimer internal energy — 
the molecule energy (1301) can be retrieved from a sim- 
pler model. A point-like boson of mass 2m interacting 
with a Fermi sea with a boson-fermion scattering length 
Oad = ao(0)a leads also to Eq. ([50)) . This means that 
the composite structure of the boson is not apparent at 
this order of the calculation. Actually this first-order 
correction stems from small values of fc ~ fc^. This orig- 
inates from the low k singular behavior in the integral 
equation and is outlined by the rescaling of variables per- 
formed during the calculation, see above. Hence, low k 
corresponds to large distances over which the composite 
structure of the boson is smeared out. 



C. Two particle- hole excitations 

The treatment becomes much more involved when 
pushed to next order in kp. First of all, Pauli block- 
ing in the two-body problem adds a new contribution to 
A. It corresponds to the ^ kp term in the expansion of 



Eq. (fTT|) . A is therefore given by 



A=l 

a 



Ix'p 
5tt 



167r- 



\v\<xf |u|>a;j!' 



where ifu^v has to be determined. There are basically 
two kinds of contribution for the solution of the inte- 
gral equation and consequently for the energy. The first 
terms come from low fc, i.e. values of k on the order of 
kp as in the last subsection. The divergence of fk,q in 



this regime implies that one has to go to higher orders 
in the calculation by including two particle-hole excita- 
tions. Nevertheless, only the k ^ part is kept which 
simplifies calculations a lot. The second kind of terms 
comes from a direct expansion of the integral equation in 
V once singular terms have been carefully treated. They 
involve larger values oi k, on the order of 1/a. 

We first need to derive a complete integral equation 
that includes all terms relevant at this order of the calcu- 
lation. The wavefunction ansatz is extended to describe 
two particle-hole excitations. It is given by Eq. ([2]) that 
we recall here. 



\ k 

k,q 



h.h' .q 



Y '^k.,k'..q..q'ai.,l(^i',lb^-k'-k+q+q'0'q,iO'q' 



(33) 



Xk,k',k", 

k.k' .k" .q.q' 



<l,<l'°^k,i°^k' .i°^k" ,i 



^^-k'-k-k"+q+q' ,]^q,i^q' A 



\\FS), 



where all q are restricted to the Fermi surface and all 
fc outside the Fermi surface. With no loss of general- 
ity, functions are antisymmetrized with respect to the 
fc variables as well as with respect to the q. This is of 
course consistent with the anticommuting properties of 
the fermionic operators. 

The Schrodinger equation is projected on the various 
subspaces of the ansatz. Eqs. (fH)) are recovered while 
Eqs. (fT6)) are changed to 



(34a) 



go + + efc H '''' j '^'' '^ 



+ 2 ^ Afc/+fc/2$fc,fc',q = — AoAfc, 

|fc'|>feF 

2 (ffc + Cfc' + e-k-k'+q ^ £q + ^k,k'.q 

+ ^k' + k/2{'fik,q - 'Pk',q) + 4:Aa Y *fc,fc',q,q' = 0, 

\q'\<kF 

(34b) 

with the notation ik = k'^/2. The k dependence of A^ 
is kept only when necessary to regularize the two-body 
problem. To get an integral equation, we proceed as be- 
fore. We obtain three new terms compared to Eq. (|25p . 
The first one is a source term proportional to ^^k.k^q.q'- 
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The second comes from the |fc'| > kp restriction in 
Eq. p4ap . Taken for fc ^ 0, it gives a term proportional 
to nF'Pk,q- The last term stems from the replacement of 
Ak using Eqs. (HH). Taken for fc ^ 0, it gives a term 
proportional to J2\q'\<kF ^k,q' ■ 

Performing the change of variable k — Xu followed by 
the rescaling ((24)) we obtain the integral equation 



1 

2^ 



\u'\>xp 



1 + + + u'u — v{u' + u) 



3u^ u ■ V 1 



47r 



4 4 2 a\ V67r2 



1 



\v'\<xp 



u'v' 



(35) 



i?* -> to 



",q,q' + '^k",k',q,q' 
A2 + fc'2 + fc"2 + k' ■ k" 



(38) 



_ fk\q' — 'Pk,q' + 'Pk,q " fk'.q 
"'^ A2 + /c'2 ■ 

The last step is to realize that the singular behavior at 
low k, (fk.q — 1/^^, implies that /c-independent terms 
become negligible for fc — > 0. Rescaling variables, e.g. 
k = Xu with X = 1/a and Eq. ((24|) , the following integral 
equation is obtained 



and the energy is given by Eq. ([32]) . 

The new corrections that we have included are all 
proportional to and our calculations goes up to 
x'jp. Therefore these corrections are negligible except for 
k kp (low k) where they are partially compensated by 
fk,q — i/kp. This justifies the fc ^ limit that is taken 
below. 

We now project the Schrodinger equation on the two 
remaining subspaces. 



1 + u'' 



(39) 



It acts only on u' and describes a three-body problem 
where u,v,v' are dummy variables. Its solution follows 
trivially from the solution of Eq. (fTO|) , 



u,u' ,v,v' 



^^•V.-^.,.')^. (40) 



ql + X"^ + Ek + Cfc' + eq - eq, I '^k,k'.q,q' 



3^ Afe//+fc//2Xfc,fc',fc",q.q' = -^($fc,fc',q' -<&fe,fc',q), 
k" 



(36a) 



3(A^ + Cfc + Cfc' + fk" + i-k-k' -k"+q+q' 



^q ^q' )Xk,k' ,k" ,q.q' 
+ Afc"+fc'/2(5'fe,fc',q,q' + ^'fc',fc",q,q' + 'I'fc",fc„q,q') — 0. 

(36b) 

These coupled equations are quite complicated. Never- 
theless many simplifications can be done that remain con- 
sistent with the order of our calculation. The q (and q') 
dependence can safely be neglected as well as the restric- 
tion outside the Fermi sea for the fc wavevectors. The 
limit fe ^ is also taken. In the source term (r.h.s.) of 
Eq. I|36ap . we can replace ^k,k'.q by using its lowest order 
expression from Eq. (j34bp . 



k,k' .q 



Aq (pk'.q - fk,. 

2 A2 + fc'2 



(37) 



Finally, Eq. (j36bp allows to express Xk,k' ,k" ,q,q' ■ Once 
incorporated into Eq. (j36ap , it leads in the universal limit 



Inserting this solution into Eq. ((35)) finally leads to the 
complete integral equation 



1 

2^ 



\u'\>X]P 



1 + u'2 + u2 + u'u — v(u' + u) 




(41) 



3ao(0) 
167r2 



d V (fu^^' = 



\v'\<Xf 



1 +U2 



Note that the l/(aA) term in Eq. ((55)) has been compen- 
sated inside the parenthesis by using Eq. ((2T)) . 

We now describe how Eq. ((41)) is solved perturbatively. 
Its solution is written in the form 



a{u, v) 



y2 _ y2 1^ _ 2u ■ V/Z 

, ao(0) 



'3„.' 



d-'v 



V ^l.'K.. u^-v'y3-2u-v'/3j 

(42) 

leading to a regularized integral equation for a{u, v) that 
can be expanded in v. The integral term inside the 
parathesis of Eq. ((1^ is a correction xp. It is im- 
portant only for small u to compensate the last term in 
the l.h.s. of Eq. (gT)). 
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Ordering the different source terms, a{u, v) can be de- 
composed as 

a{u, v) = asi{u, v) + as2{u, v) + ai{u){u ■ v) + u^a2(u), 

(43) 

wliere tfie different terms are to be detailed below. 
asi{u, v) is the same as in the last subsection and is given 
by Eq. (|^^ . Its lowest order in xp is ao(u), solution of 
Eq. (fT9|) . The first correction in xp has been calculated 
last subsection. The expansion of Eq. ([^^1) also yields a 
Xp correction. Once plugged into Eq. ([32]) . aX receives 
the correction —aQ{0)xpl3 with 



3 

47^3 



[ dss^ [2 -sJi(s)]2 -0.01419, (44) 
Jq 

where Ji(s) is given Eq. (|28p . The second term, (u), 
corresponds to the last term in the l.h.s. of Eq. (|¥T|) . 
Similarly to asi(u), it stems from low k kp, and re- 
duces, up to an important prefactor in the source term, 
to the three-body problem integral equation The 
solution can be written as 

2 



as2{u,v) = - 
where 

J2{v) = / d'x- 



Qo(0)a:F 
27r2 



aa{u)J2 ( — ) , (45) 



1 



x>l X" 



■2 o , d^v' - ,1 o , • 

y£ _ 2x y / , ™2 _ j/1 _ 2a; y' 

3 3 Ji/ <i J- 3 3 

(46) 

The contribution to a\ is given by a^{^)x\l2 with the 
integral 



1 



167r5 



dyy^J2{y) 0.014440. 



(47) 



In a way similar to the first order correction, the two 
contributions stemming from a^i and as2 are produced 
by small wavevectors k kp and match exactly the total 
x'jp correction of the point-like boson model. They cor- 
respond to the large distance part for the scattering of 
the boson by the Fermi sea. On this scale k ~ kp, the 
composite structure of the boson is not apparent. 

This is in contrast with the situation for the last two 
terms in Eq. They result from the direct ex- 

pansion in V of the integral equation. In that case, 
typical wavevectors k are of order 1/a 3> kp. These 
terms involve the whole spatial structure of the composite 
molecule. They have no equivalent in the point-like bo- 
son model and hence reveal the composite nature of the 
molecule. The integral equations that determine ai{u) 
and a2{u) follow from a tedious but systematic expan- 
sion in V. 

^+00 

du'Ki {u,u')ai{u') = Si {u) (48) 



determines ai{u) with the Kernel 
i5{u-u') , lu'^ 

Ki(u,u) - 



4 l + E 









hi- 











(49) 



with A = l + v? + B = uu' and E = ^1 + iu^/A. 
Using the same definitions, the source term can be writ- 
ten 



Si{u) 



aoiu) 



" 



TT 



4ln 





1)] 







The integral equation for 02(1*), 





^A + 




hi- 











(50) 



+ 00 



du'Kn{u,u')a2{u') = s!^^\u) + S^-"' (u), (51) 



(2), 



has the same Kernel as Eq. ([T9| with two sources terms. 
Note that this result has already been averaged over the 
directions of v and u since only that component matters 
for A, see Eq. ([5^ . 

The first source term, S2^\u), depends on ai(u) 

(1) _ 1 u^aiju) 
^2 16£;(l + i?)2 



2 

3^ 
2 

97r 



+00 



du'ai{u') 



+00 2 
du' ai(u') — In 



■In 



A + B 
A-B 



(52) 



The second source term depends on ao(u), 

3ao(u) u'^{l + 3E)ao{u) 



5f N = - 



32^(1 + .B)2 128E3{1 + E)^ 



du' 
12 



aoiu') ^^f A + B 



.A-B 
P-A 



26 

27^ 7o "'^ L 25 

4 /•+°° ,a^{u') v!--A ^ 1 

2 , (u^W^)^^ 



1+U2 





^)1 







(53) 



The resulting contribution to oA, Eq. (|32|) . is added to 
the Pauli blocking term (third term inside the parenthesis 
of Eq. ([5^ ) to give with 



h 



3a2(0) 3ao(0) 



57r 207r 

g /-Hoc 



3207r 157r2./o (1 + ^2)2 







flo(M)(l + 5mV2) 

(1-|-U2)2 



0.063324. 



(54) 
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Gathering all xj^ corrections to the mean field result, 
we finally obtain that Eq. ([50]) is extended to 



E = -- 



l-%(^^^->|W(ln2-3/8) 



+ 2a;^[a^(0)(/2-/3)+/i] , 



(55) 



with Eqs. ([Ml) , (SZl) > iMl) for h^h and h- Note that 
ao(0)(/2 - h) + h ~ 0.0637. It indicates that the con- 
tribution Ii, specific to the composite boson, dominates 
quantitatively the last term in the energy. The point-like 
boson model also leads to Eq. ([SS]) but with /i = 0. The 
final result ([55)1 was announced in Sec.|lT]as Eq. 



IV. VARIATIONAL TREATMENT 

As mentioned in Sec. HIl where the main results of this 
article are summarized, an alternative way of treating the 
problem is to use a variational approach in the restricted 
subspace given by Eq. ([4]) . We recall here the variational 
form 



t 

-fc,T 



(56) 



where ^k,k' has been antisymmetrized. As observed in 
the previous section IIIIl this ansatz is sufficient to re- 
cover the mean-field correction pip corresponding to the 
first order term in the perturbative expansion (j55p . Min- 
imizing the energy with respect to the amplitudes /3, Ak , 
ipk and $fc,fc', while keeping constant leads to a 

set of equations very similar to the one obtained in the 
perturbative treatment of the problem 

(Eo - E)(3 + -y^KkAk = (57) 

k 

{2ek - E)Ak + + -^Y^^Wk = (58) 

V V q 

Xl'(^o + Efe + eq-fc/2 -Eg- E)ipk + 



q 

A 



q k' ,q 

E ('^'^ + '^k' + £q-k-k' -£q- E)^k^k' + 

q' 

^y^E (^fc'+fe/2 Vk - Afc+fe72 Vk') = 



(59) 



(60) 



Note that here the energy E appears as the Lagrange 
multiplier associated with the constraint on {iplip). More- 
over, since the equations are linear, we can set /3 = 1. 

The first two equations are identical to those we used 
in the perturbative calculation. Introducing A defined by 
E ~ —h^X^/m, it yields the first variational equation: 



^ _ 1 2kf 
a TT 



1 — — arctan(A:i;-/A) 



k 



(61) 

In Eq. (|60|) . the explicit q dependence vanishes with 
the radial sum, and we finally recover the equation 



2^V 



^k + £k' 



£k+k' — E 



Inserting Eq. ([55)) and (p^ in ([55)1 finally yields in the 
short range limit a closed equation for the field (p: 



where the function F and the integral kernel L' are de- 
fined by 



Vk' 



1 



(64) 



'V^ E 
q 



2£k 



1 

47r 



A2 



1 

V 



E 

|fc'|<fcF 



£k + £fc' 



£fc+fc' 



E 



(65) 



Eq. (|6ip and (|63p are then solved numerically, yielding 
the result presented in Fig. ([3]). A good agreement be- 
tween our variational ansatz and the exact Monte-Carlo 
simulations is obtained. As already noted, it is very close 
to the mean-field prediction up to kpa ^ 1. However, it 
stays finite at unitarity, since we have E ~ 1.1506 for 
kpa = oo. 



V. CONCLUSION 

In analogy with the Fermi-polaron system, the work 
presented here shows that the molecular sector of the 
impurity problem can be described quantitatively as a 
molecule dressed by a single particle-hole excitation. 
However, two important points are still to be clarified 
and will be addressed in future work. First, what is the 
effective mass of the quasi-particle ? This property is in 
particular important to capture the dynamical behavior 
of the system, as suggested [l4| and observed recently in 
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experiments [23j. Moreover, the bosonic/fermionic na- 
ture of the quasi-particle should be clarified by the study 
of an ensemble of impurities immersed in a Fermi sea, 
that may help clarifying the molecule-polaron transition 
scenario. This formalism can also be used to interpret 
spectroscopic data obtained for instance in Ref. [131 and 
refine the Nozieres Schmitt-Rink analysis of Ref. [25| . 
Another extension of this work deals with mixtures of 



particles with unequal masses, for instance the Li-K mix- 
ture. 
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